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ABSTRACT 

A nearby star having a near-transit of a galaxy will cause a time-dependent weak 
lensing of the galaxy. Because the effect is small, we refer to this as weak microlensing. 
This could provide a useful method to weigh low-mass stars and brown dwarfs. We 
examine the feasibility of measuring masses in this way and we find that a star causes 
measurable weak microlensing in a galaxy even at 10 Einstein radii away. Of order 
one magnitude I ^ 25 galaxy comes close enough to one or other of the ~ 100 nearest 
stars per year. 
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1 INTRODUCTION 

Gravitational lensing by Galactic stars has come a long way 
since the low-probability assessment of |Einstein| (|1936| ) . The 



first detections (Alcock et al. 



1993 



Aubourg et al. 



1993 



Udalski et al. 1993) have been followed by ~ 4000 more, 



including a few detections of planets around the main lensing 
star (e.g.,|Gould et al.||2006|). 



Yet despite the now abundant examples of microlensing, 
accurate measurements of the lensing mass are still rare. The 
reason is that a microlensing light-curve on its own tells 
us precisely how long proper motion takes to traverse an 
Einstein radius, but the physical length te of the Einstein 
radius and its angular size 0e remain unknown. Hence, only 
a broad statistical statement about the stellar mass can be 
made in most cases. 



Refsdal (1966) anticipated the problem and suggested 



a way to overcome it: If both stars are directly observable, 
their relative proper motion would (together with the light 
curve) supply 0e, while a second observatory in the solar 
system would provide te- Several related ideas have been 
forwarded, notably by |An et al.| ([2002 ), who exploited a 
combination of caustic-crossing times and finite source size 
effects to obtain the first microlensing mass measurement. 
A similar strategy has recently been used by | Gould et al.| 
(2009) to measure the mass of a brown dwarf to 10%. 

Another possibility is to measure the lensing effect of 



a star on a galaxy. Paczynski ( 1996 ) pointed out that as 



a Galactic star makes a near-transit of a galaxy, the lat- 
ter will undergo a small shift in its apparent centroid. For 
this variety of microlensing to be observable, 0e needs to 
be larger than the precision of image centroiding-although 
it can be below resolution-and proper motions need to be 
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large. Nearby stars are the only realistic prospect, since 
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(1) 



for sources at infinity, and proper motions are ~ \" per year. 
Microlensing by nearby stars would have none of the degen- 
eracy problems mentioned above; the lens being at known 
distance, the mass is the only unknown parameter. More- 
over, such events are predictable long in advance. The transit 
may be on the order of a few months allowing for observa- 
tions of the galaxy both before and during the lensing event. 
Ideally, the star would not transit directly across the galaxy 
since the star must be masked out for proper observations 
of the galaxy. 



Paczynski ( 1996 ) suggested that microlensing centroid 



shifts could be used to measure down to masses of nearby 
brown dwarfs. In this paper, we suggest a refinement of 
Paczynski's idea which could make it much more effective. 
Rather than just the galaxy centroid shift, the whole weakly- 
lensed image of the galaxy could be exploited to infer the 
lensing mass. We develop a technique to extract the weak- 
lensing effects and estimate the mass of the star. 



2 A FITTING METHOD 

Consider a star with Einstein radius Be whose sky position 
at epoch t is zt. Lensing by this star maps a point 6 in the 
image plane to a point <j) in the source plane such that 



(2) 



where zt,0,<f> are two-dimensional vectors. 

Next, we consider a galaxy whose unlensed brightness 
distribution S is expanded in terms of basis functions B n (0) 
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S{6) = ^2a n B n {6) 



(3) 



where a n are the expansion coefficients. In the presence of 
lensing, the brightness distribution of the galaxy will be 



s t {e) = Y,a n B n (<t>(6,z t ,e E )) 



(4) 



where S now has a time dependence due to the position of 
the star zt- Since lensing conserves surface brightness, the 
lensed surface brightness at equals the unlensed surface 
brightness at (j>. 

Pixelating the image plane, we write the pixel- wise 
brightness distribution as 



(5) 



where L™ ^ represents a lensed and pixelated basis function 



(6) 



The expression in Equation [5] is our model for the lensed 
brightness distribution. If the observed distribution is dt,ij 
then the likelihood is 

C{a n ,0 E ) oc ]^[exp [-§ cj~ 2 %3 (d t ,ij - D t ,ij) 2 } (7) 

where a t ,ij is the pixel- wise noise, which we assume is 
Gaussian. The model has a complicated dependence on fe, 
but only a linear dependence on the expansion coefficients 
a n . Since we are not especially interested in the a n , we 
marginalise them out by standard methods (for example, 
chapter 5 of Saha 2003). The marginalised likelihood is given 
by 

2 In C{6 E ) = In | det C\ + ^ PmP n C mn - °t,l d la • 



t,ij 



Here 



Pn = &t,ij dt,ij L™i 



(8) 



(9) 



t,ij 



represents a kind of projection of the data on the model, 
while 



— 1 \ ^ —2 j m jn 

mn = / J a t,ij ^t,ij ^t,i, 

t,ij 



(10) 



is the inverse covariance matrix. We can relate C(6e) to an 
effective \ 2 as j us t 



C(0 E ) = ex P (-x 2 /2) 



(11) 



We are now prepared to estimate the mass of a star by 
its lensing effect on a background galaxy. Given the pixel- 
wise brightness dt,ij and noise level cr t ,ij, we simply need to 
calculate the effective x 2 and minimize with respect to 0e- 

As an aside, | det C| in Equation [8] will typically over- 
flow floating-point arithmetic, while ln|detC| will fit quite 
nicely. Hence we compute an alternative form, namely, 
In | det C\ = m ^n, where A n is the nth eigenvalue of C. 



7tot |,1,3,7,12 xlO 6 

A^obs 2,3,4,5 

9 E , true 10, 20, 30, 40, 50 [milliarcsec] 

p 0.08, 0.126, 0.232, 0.454, 0.903 [arcsec] 



Table 1. Parameters used for simulated observations: jtot is the 
total number of photons over all observation epochs (including 
one unlensed observation), A^ Q b s is the number of epochs, #£;,true 
is the Einstein radius we wish to recover, and p is the closest 
approach of the star on the sky plane to the background galaxy. 



3 SIMULATED OBSERVATIONS 

We tested the above scheme using simulated data, focusing 
specifically on the dependence on four quantities: the actual 
Einstein radius of the star e, true, the closest approach or 
impact parameter p, the number of observed epochs (or t 
values) iVobs, and the total number of photons collected 7 to t- 
The values are summarized in Table [T] The full matrix of 
these parameters was tested, giving a total of 500 simulated 
observational programs. 

For the exact form of the unlensed surface brightness in 
Equation [3] we chose 

S(0) = exp (-7.67 yie/ReY + R^ V * j (12) 

with R e — 1. This is simply a de Vaucouleurs profile mod- 
ified by a core radius R c of 2 pixels to mimic the effect of 
the telescope PSF on a singular cusp. 

We then considered 71 x 71 pixels imaging a patch of 
sky 2" on a side and centred on the galaxy. With this reso- 
lution each pixel is about / = / . / 028 across, which, for exam- 
ple, would be equivalent to about one pixel per resolution 
element of the Nasmyth Adaptive Optics System (NAOS) 
Near-Infrared Imager and Spectrograph (CONICA) camera 



(Lenzen et al. 2003 Rousset et al. 2003). At a redshift of 



z ~ 0.5, where we expect to find most of our background 
galaxies, \" is about 6 kpc. The corresponding lensed pixe- 
lated brightness is 

(13) 



St,ij — S((j)(0ij, Zt, 9 E, true)) 



We normalized this brightness to the total number of pho- 
tons 



&t,ij — 7tot, 



(14) 



and thus 7tot/Aobs is the number of photons per image. Con- 
sidering the normalized St,ij as the expectation value of the 
pixel- wise photon count, we then drew the simulated data 
from a Poisson distribution. The noise level at.ij was 



H,ij 



i — rr n 

taken as ydt,ij, which greatly simplifies Equations M and ^1 
In practice, the star will be masked out during observa- 
tions. To simulate this, we discarded the pixels within one 
pixel- width of the star. 

As basis functions we chose two-dimensional Hermite 
functions or shapelets (e.g., Refregier 2003). The scale pa- 
rameter for the shapelets was set to f3 = 0^2 and we used 
Abasis = 20 x 20 basis functions in all. These settings allow 
structures as large as # ma x = /3(2\/Nb&sis + 1) 1//2 = 1//28 
and as s mall as 6> min p(2^ N hasis + 1)" 1/2 = 0"03 to be 
resolved ( |Melchior et al.|2007 ). 

For each of the 500 test cases, we considered one epoch 



Weak Microlensing 3 




Original 



Reconstructed 



Residual 



0.05 0.06 
9e [arcsec] 

Figure 1. Plot of the effective x 2 as a function of Oe, for a sim- 
ulated observation program with N \> s = 2 epochs, total photons 
7tot = 3 x 10 6 collected, #£,true = 0'/05 and impact parameter 
p = / /126. The solid vertical line marks e, true and the dashed 
line marks the best fit. The formal count of degrees of freedom is 
2 epochs x 71 2 pixels — 20 2 basis functions = 9682, so the effective 
reduced % 2 ~ 1.4. The dependence of the fit on Be is non-linear, 
hence the asymmetric shape of the curve. 



t = with the galaxy unlensed and additional epochs t > 
with the star at z t = (bit, -12* \ i e [1, . . . , 5]). The selec- 
tion of i corresponds to the selection of p: A large p implies 
a large i. This choice of position puts the star at the centre 
of a pixel. Extensive testing has shown a sensitivity to plac- 
ing stars near pixel edges, whereby recovery of the data is 
degraded if the star is too close to pixel boundaries. 

In Figure [l] we show the effective % 2 as a function of 0e 
for one of the simulated data sets. In Figure [2] we show the 
simulated images, along with the reconstructed and residual 
images for the best-fit Be- Examining such plots is a good 
indicator of potential problems. For example, if too few ba- 
sis functions are used, a grid- like pattern shows up in the 
reconstructions and the residual, and recovery of 0e is very 
poor. 

Figure [3] summarizes the complete suite of 500 simu- 
lated observation programs, showing the mass-recovery er- 
rors as a function of #E,true, P, 7tot and iV bs- The following 
conclusions can be easily read off: 

• The mass range of nearby brown dwarfs is accessible, 
since E down to 0"02 can be measured with the resolution 
considered. 

• Impact parameters of p ~ 10 Be are small enough, but 
if p is too small the galaxy can be obscured by the star mask, 
leading to poor results. 

• Of order a million photons from the galaxy are needed, 
and a few times this are desirable, but it does not matter 
much whether these are concentrated in two epochs or dis- 
tributed among several epochs. 



4 EVENT RATES 

We now consider how likely is it to find a star crossing near 
a background galaxy. For this analysis we used the Research 




Figure 2. Details of the best-fit 0e for the simulated observation 
program of Figure ^ The upper row refers to an unlensed epoch 
t = and the lower row to the lensed epoch t = 1. The left column 



shows the simulated data d, 



the middle column shows the best 



reconstruction Dt,ij, while the right column is \dt,ij —D t ,ij\. The 
scale is in units of photons. The star has been masked out as 
would be done during an observation. 
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Figure 3. Mass uncertainty as a function of four quantities: 0e, 
impact parameter total galaxy photons collected 7tot , and the 
number of epochs Af b s - Within each box, 0e and p are varied 
at fixed 7tot,^V b s - The latter two quantities are varied between 
boxes, as labelled. Circles indicate the error in mass (i.e., in E )\ 
Large, green circles denote < 5% error, orange circles 5-11% error, 
red circles 11-20% error, and small, open circles are used if the 
mass error was > 20%. Missing circles mean that no mass estimate 
could be made. Of all the tests, 26% have errors less than 5% and 
39% have errors less than 11%. Of those tests with filled circles, 
52% have errors less than 5% and 79% have errors less than 11%. 



Consortium on Near by Stars (RE CONS) list of the 100 near- 
est stellar systems ( Henry |[2009 ) . The proper motions and 
estimated masses of the stars in these systems are plotted 
in Figure [4] In Figure [5] we have plotted the area of sky 
swept out by Einstein radii per year, or 2 \iQe where \i is the 
proper motion. 

If we restrict ourselves to masses < O.5M and proper 
motions > ; /5 /yr, we are left with 85 stars. Assuming, as 
seen in our tests, that a galaxy within 10 Be is a candidate, 
we sum 20 /jlOe over these stars. The total available area is 
per year. 



The GalaxyCount program (Ellis & Bland-Hawthorn 
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Figure 4. The nearest 100 stellar systems (139 stars) from the 
RECONS catalogue. Shown are the proper motions (left) and 
estimated masses (right) for the sample. Systems in red have been 
excluded from the discussion in Section^] The five highest mass 
stars are, a Centauri A+B, Sirius, Procyon, and Altair. 



brown dwarfs with a mass of ~ 0.05M© (0e ~ Of! 01) such as 
SCR 1845-6357 B (at 3.9 pc), DENIS 0255-4700 (5.0 pc), 
2MASS 0415-0935 (5.7 pc), or GJ 570 D (5.9 pc), have 
I ^ 17 — 20 mag and V ~ 22 — 25 mag and they are not or 
not much brighter than the abundant backgound galaxies. 
Low mass stars and substellar objects are red or extremely 
red and imaging observations of blue star-forming galaxy at 
short wavelengths is favoured because the image contamina- 
tion of the lensed galaxy by the PSF of the lensing object 
is strongly reduced. It seems that HST or an adaptive op- 
tics systems (e.g. with laser guide star) at a large telescope 
working at short wavelengths < 1/xm should be capable of 
achieving successful observations for certain weak microlens- 
ing events. 
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Figure 5. Sky area coming within one Einstein radius of each 
nearby star per year, plotted against distance (left), and against 
estimated stellar mass (right). Systems in red have been excluded 
from the discussion in Section 0] 



2007) estimates c± 1 galaxy with magnitude I ^ 25 within 
a sky area of 70arcsec 2 . This provides a rough estimate of 
the rate of observable weak microlensing events. 



5 OBSERVATIONAL PROSPECTS 

Observing the weak lensing of a faint galaxy by a nearby star 
would require a high resolution (« Of! 05) imager with high 
contrast capabilities. A O.5M star at 5 pc has a brightness 
of I ~ 6.5mag while I ~ 12mag for a 0.1 Mq star at the same 
distance. Thus a contrast in the range Al — 12 — 18mag must 
be achieved at a separation of about 0"05 for typical events. 
This is quite a challenge for existing instruments. Fortu- 
nately, rapid progress can be expected in this field by instru- 
ments currently built for the imaging of planetary systems 
with 8m- 10m telescopes and further significant progress will 
be possible with extremely large telescopes and high con- 
trast imagers in space. They will provide very high contrast 
observations Am > 20mag and allow mass determinations 
of many nearby stars using weak microlensing of faint back- 
ground galaxies as advocated in this Letter. Any background 
light that is not from the galaxy can still be considered part 
of the source as it will either be lensed or remain relatively 
constant throughout the duration of the complete observa- 
tion program. An 8m class telescope with 30-50% efficiency 
collects about 50,000 photons/hr. Thus, a typical program 
might need between 20 and 100 hours to expect reasonable 
results. 

With existing instruments it should already be possi- 
ble to observe weak microlensing in favourable cases where 
the impact parameter is small and the optical resolution 
is higher than that considered here. Nearby (d ~ 5 pc) 



REFERENCES 

Alcock C, Akerlof C. W., Allsman R. A., et al. 1993, Na- 
ture, 365, 621 

An J. H., Albrow M. D., Beaulieu J.-R, et al. 2002, ApJ, 
572, 521 

Aubourg E., Bareyre P., Brehin S., et al. 1993, Nature, 365, 
623 

Einstein A., 1936, Science, 84, 506 

Ellis S. C, Bland-Hawthorn J., 2007, MNRAS, 377, 815 
Gould A., Udalski A., An D., et al. 2006, ApJ, 644, L37 
Gould A., Udalski A., Monard B., et al 2009, ApJ, 698, 

L147 

Henry T. J., , 2009, RECONS, |http : // j oy . chara . gsuT] 

|edu/REC0NS71 

Lenzen R., Hartung M., Brandner W., et al. 2003, in 
M. Iye & A. F. M. Moorwood ed., Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Se- 
ries Vol. 4841 of Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, NAOS-CONICA 
first on sky results in a variety of observing modes, pp 
944-952 

Melchior P., Meneghetti M., Bartelmann M., 2007, A&A, 
463, 1215 

Paczynski B., 1996, Acta Astronomica, 46, 291 
Refregier A., 2003, MNRAS, 338, 35 
Refsdal S., 1966, MNRAS, 134, 315 

Rousset G., Lacombe F., Puget P., et al. 2003, in 
P. L. Wizinowich & D. Bonaccini ed., Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Se- 
ries Vol. 4839 of Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, NAOS, the first AO 
system of the VLT: on-sky performance, pp 140-149 

Saha P., 2003, Principles of Data Analysis. Great Malvern, 
UK, Cappella Archive, 2003 

Udalski A., Szymanski M., Kaluzny J., Kubiak M., 
Krzeminski W., Mateo M., Preston G. W., Paczynski B., 
1993, Acta Astronomica, 43, 289 



